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We show that the entanglement entropy associated to a region grows faster than the area of its 
boundary surface. This is done by proving a special case of a conjecture due to Widom that 
yields a surprisingly simple expression for the leading behaviour of the entanglement entropy. 
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1 Background 

An interesting model for the area law of Bckcnstein's black-hole entropy is a local quantum field 
theory for which one restricts all observations to the complement of a compact spatial region fi. 
One should think of f2 as the interior of the horizon although one does not require O to have any 
gravitational relevance. In fact, we will restrict our attention to a region in the n-dimensional fiat 
Euclidean space R". Even if the quantum field is globally in a pure state (e.g., the vacuum), the 
• state restricted to the complement of Q will in general be mixed with finite von Neumann entropy, 

| called the entanglement entropy, S(£l) [TJ[2J[3]. For many types of quantum field theories, it was 

f~>) ■ observed that in the semi-classical limit of large fl, the entanglement entropy S(fl) scales like the 

f^*) | area of dtl. In other words, if we rescale fl C M™ by a factor R, the entanglement entropy, S(RQ), 

should asymptotically scale like i?' 1 " 1 for large R. 

Recently, it was observed by Gioev and Klich HJ|S], that for free fermions at zero temperature, 
however, the entanglement entropy of the ground state scales like log(i?) for R — > oo based 

rS | on a conjecture by Widom [6l [7]. They provide a lower bound on the entanglement entropy in 

terms of the trace of a quadratic function of the restricted state. For the latter, we have proved 
in [8] the leading asymptotic behaviour. This establishes indeed that the entanglement entropy 
scales at least like i? n_1 log(i?), which violates of the beforementioned area law scaling. 

For notational clarity we will often write equalities that only hold asymptotically for large R, 
that is, we will drop sublcading terms that are not central to our argument and only mention it 
in the text. 

General arguments imply that for a pure state, the entanglement entropy of this state with 
respect to fl is the same as for the complement M. n \ fl. For simplicity, we will restrict the state 
to a compact region O. The ground state, pr, of a system of non-interacting fermions in W 1 is 
given in terms of the Fermi surface at Fermi energy ep'- All one-particle states with momentum 
p G T = {p e R n |i?(p) < cf} and energjQs(p) are occupied. The ground state pr is then 
characterised by the one-particle Fermi projector, Pr, defined in momentum space by the kernel 
-Pr(PjP') = Xr(p)<5(p — p')- Here, and in the following, \A denotes the indicator function of a 



O 

m 



1 We will not discuss here assumptions on the dispersion relation -B(p) but one may, of course, think of the 
example -B(p) 



Copyright line will be provided by the publisher 



2 



RC Helling and W Spitzer: Fermions Violate Area Law 



set A. In position space, the kernel Pr(x, x') = x?( x — x') with xr(x) = (27r) _Tl J r <ipe lx ' p being 
the inverse Fourier transform of xr- 

In order to restrict the state pr to the region St wc project the Fermi projector Pp onto L 2 (Sl) 
with — Xfi- This gives the reduced one-particle density matrix £?n,r = QnPrQn- The 
entanglement entropy, S(f2, T), of the many particle system in the ground state pr restricted to 
the region Si is then defined as the grand canonical entropy of £>n,r, that is, S(fl, F) = tr 77(^0^) 
with rj(t) = -Hog (t) - (1 - 1) log (1 - t) for < t < 1. For details see Section 4]. We are 
interested here in the behaviour of this entropy for fixed T but large SI. To this end, we also fix 
SI and study the asymptotic behaviour of S(RSl,T) as R — > 00. 

Our main result is the asymptotic computation of tr[pfln,r(l — £>rai,r)] as R — > 00: 

tr[ ejai ,r(l - Qan,r)] = (^)" ' f dA(x)dA(p) |n x • n p |, (1) 

up to terms that grow slower in R. Here, n x denotes the unit normal vector at x £ dSl, cL4(x) is 
the surface measure on 917, and similarly for n p and dA(p). 

Since r)(t) > log(2)4i(l — t) for < t < 1, we obtain the asymptotic lower bound on the 
entanglement entropy, 

S(M1,T) > log(R) f dA(x)dA(p) |n x • n p |. (2) 

Gioev and Klich conjectured in [5] that the exact scaling of S(R£l, T) is obtained if we replace the 
factor log(2)/7r 2 in ([2]) by y^. This remains an open problem. 



2 Computation of tr[^n,r(l - Qm,r)] 

The trace of Qrqx is simply equal to (^:) n |f2||r|, where | ■ | denotes the n-dimensional Lebesgue 
volume. The trace of g Rn r equals 

tr(Q R nPrQRnPr) = [ d* [ dad |x?(x - x')| 2 = / dv |x?(v)| 2 \RQ n {RSI - v)| . (3) 
Jrq Jrq JR(n~n) 

In the last line we have changed the variables x and x' to u = x and v = x — x' . Then we expand 
the volume \RQ n (RQ - v)| to first order in v (cf. [9j Theorem 2.1]), 

\RSID (RSI - v)| = R n \n\ +R n ~ 1 [ dA(x) max(0, v • n x ) + i?"- 2 0(|v| 2 ). (4) 

Jan 

Let us first look at the contribution of i?"|f2| to the trace of The function v i-> x?(v) 

decays like |v| _( ^™ +1 ' / ' 2 for large |v|, see ([7]). At the cost of an order i? Tt_1 -term we may therefore 
extend the v- integration to all of R". By the Plancherel formula this integral gives (27r) _ ™|r| and 
cancels with ti qrut- When integrated over v, the remainder term R n ~ 2 0(\v\ 2 ) is also seen to 
yield a term of the order R n ~ l by using again the above mentioned decay of xf. 

Thus, our bound on the entanglement entropy will come from the second term in ((4]). Here, we 
write max(0, v ■ n x ) = X[o,oo) ( v ' n x) v ■ n x and use the Gaufi Theorem so that 

(27T)" v x?(v) = -i [ dA(p) n p e lv 'P. (5) 
J or 

So it remains to show that for p G dT, 

rfvx [ o,oo)(v-n x )xB(-v)e lv -P + (27r l )- 1 sgn(n x -n p )log(i 1 '.) = o(R) . (6) 



'ij(n-n) 
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Let us now consider the function v i— > xr(— v) in detail. We use the representation from ([5]), 
that is, (27r) _n xr(— v) = ■ J gr dA(p') p' e _tv ' p . Then we introduce a coordinate system 
where v = (0, . . . , 0, v) and where the boundary <9L is locally written as the graph of a function 
/: U C R n_1 -> K, that is, p' = (t, /(t)) and dA(p') = y/l + |V/| 2 dt. The unit normal vector is 
n P < = sgn(v • p')(-V/, 1)/V1 + |V/|. Then, 

/ dA(p') (rv)n e- ?v ' p = -- / ^ sgn(/(t)) e~ iv ^\ 
Jf(U) v Ju 

In order to find the leading asymptotic behaviour of this t-integral for large v we apply the method 
of stationary phase. Let k a — k a (v) = (t a (v), /(t a (v))) be the collection of all stationary points 
of such local functions /, that is, V/(t a ) = 0; in other words, the points k a S dT are such that 
the unit normal vector rik a at k a is parallel to v. Thus, 

xB(-v) = -*(2t™)-(" +1 >/ 2 Y ' gn(v • ka) e-iv^-ii ■«»(/«(*-)) + ( w -(«+D/ 2 ). (7) 

f l3 {t a )\ 

Here, fij(t a ) denotes the Hessian of / at t a , sgn(/y (t a )) the signum of this Hessian. The deter- 
minant, det/jj(t a ), equals the Gaufiian curvature of dT at k a . 

Using ([7]) we return to the oscillatory integral in ©, and employ once more the method of 
stationary phase. The composite phase from ([5]) and ([7]) is equal to v • (p — k a (v)). Next, 
we introduce generalised spherical coordinates for v as v = p(u,h(u)), where the map h: V C 
R n_1 — > R locally parametrises the boundary d(fl — f2) 3 (u, h(u)) and p £ [0, R] is a radial 
coordinate. The stationary point k a is now a function of u. We have the freedom to assume that 
n p = (0, . . . , 0, 1) such that k a ((0, h(0))) = p for one index a and that h has its extremum at the 
origin u = 0. Instead of integrating p from to R we will actually integrate p only over [1,R] 
thereby making an irrelevant error independent of R. 

To express k a (v) as a function of u we equate v/v = n ka ^ v - ) , that is, 

=sgn(v.n k J ^ . (8) 

v^+W Vi + iv t /(t)| 2 

Taking derivatives and evaluating at u = 0, we find g^-(O) = — /^ 1 (k a (0))//i(0), where k a (0) is 
short for k a (0, /i(0)) and /j" 1 is the matrix inverse of the Hessian of /. With this, we can expand 
the phase of the remaining v-integral to second order as 

f~ 1 (k (0)) 

v • (p k (v)) = ph(0)(p k Q (0))„ + P ,J 2fe( ° u iUj , (9) 

The volume element is given by dv = p n ~ 1 h(x\)dpd\i. The stationary phase integral over the 
n - 1 coordinates u yields {2h{0) / p)^-^' 2 yj\ dot fr. 1 1 exp(iph(0)(p - k o (0))„) + i\ Bgn(/y(t a )). 

This, surprisingly, reduces the remaining p-integral to dp e vfi ( " p ^ ka ( )'"//). As noted above, 
for one index a, the exponent vanishes as p = k a (0) and the integral is the desired log(i?). If, 
however, p ^ k a (0) the integral is bounded for large R and does not contribute to leading order. 

Collecting all terms we have thus proved the lower bound ([1} which grows faster than the area 
law scaling i? n_1 by a factor of log(i?). Note that the stationary phase integrals localises v and 
p' such that there are only contributions from p, n v , n p > all being parallel. 



3 Discussion 

Above, we presented a lower bound for the entanglement entropy which violates an i?™ -1 area 
law scaling. The expression ^ is in terms of integrals over the boundaries dfl and dT and is thus 
still a "boundary effect" . 
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We find it most curious that although we have used stationary phase methods, which in general 
depend on (the existence of) second order derivatives at the stationary points, all these curvature 
terms involving and hij eventually cancel out and the integrand in the bound does not contain 
derivatives. However, one expects [H [10] that a fractal boundary dVl of dimension n — 1 + a with 
< a < 1 leads to a scaling of the entanglement entropy of at least of the order R n - 1 + a and, 
presumably, without a log (R) correction. 

For n = 1, if Q, and T are a disjoint union of k and £ compact intervals of finite length 
respectively, then our lower bound (J2j) for the entanglement entropy gives log(U) log(2)4H/7r 2 . In 
this one-dimensional case, the precise scaling has been proved, namely, log(R)k£/3. If we then 
consider the hypercubes, say = T = [0, 1]™, it is not difficult to derive the exact asymptotic 
scaling of S(RCl, T) to be (i?/(27r))"^ 1 log (R)n 2 /3, which is in agreement with the conjecture by 
Gioev and Klich mentioned at the end of Section 1. Our method of proof requires that the surfaces 
d£l and dT are C 3 . Hence, hypercubes are not included. For a C 3 surface, 9r, it was crucial that 
Xr (v) behaves like |v|~(™ +1 )/ 2 for large |v|. This is not the case for a non-smooth surface such as 
the hypercube, where X\p^(v) = 117=1 sin(vi) / (nvi) . 

It should be noted as well that the discontinuity of xr is crucial for the decay of its Fourier 
transform. For example, for the equilibrium state at positive temperature, the entanglement 
entropy scales like the volume R n , see [5]. We did not have to introduce an ultraviolet regulator 
since momentum integrations are limited to the compact region T. 
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